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Abstract 


A general kinetic framework to study epoxy polymerization is 
developed. Best-fit values of five rate constants are obtained using some experimental 
chromatographic data available in the open literature on a few low molecular weight 
molecular species. Detailed sensitivity studies are then carried out to identify the most 
important rate constants. Average molecular weights and the polydispersity index are 
predicted using these parameters. The present work is more general than earlier 
kinetic models, and does not have the drawbacks of the probabilistic models available 
in the open literature. It can easily be extended to model industrial reactors. 
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Chapter 1 
Introduction 


Epoxy polymers are formed by the reaction of an epoxy group f-CHCHi) with 


a hydroxyl group. These resins are manufactured using two common processes, the 
taffy process and the advancement process, but adaptations of these are not 
uncommon. In the taffy process, a diol, such as bisphenol A lHO~ R-OH, 
where R is Q,Ht C(CHh)f 0 ,Hs 1 is reacted with a controlled excess of 
epichlorohydrin (^Tf^HCb^CI) give a mixture of low molecular weight 
O 

oligomers having glycidyl ether groups (- OCH 2 <^H^H 2 ) at both ends. The lowest 

O 

molecular weight oligomeric product is referred to as the diglycidyl ether of 
bisphenol A (DGEBPA). This process is carried out under the catalytic influence of an 
alkali such as sodium hydroxide. Water is present in the heterogeneous reaction mass, 
and it is claimed 1 that the reaction takes place at the water/organic-medium interface. 
In the advancement process, bisphenol A is reacted with equimolar amounts of 
commercially available DGEBPA. This is a typical AA + BB step growth 
polymerization, 2,3 in which complications are introduced by the presence of (over 10% 
of ) impurities, e.g., higher molecular weight diepoxides, monofunctional epoxides, 
etc. 4 ’ 5 Yet another process 5 involves the addition of bisphenol A to a solution of 
epichlorohydrin in a solvent (e.g., methyl cellosolve). Pulverized sodium hydroxide is 
added to this reaction mass in a controlled manner over a period of time. A specified 



temperature history is effected and intermediate washing with water and distillation 
under reduced pressure are used to give the final resin. 

It is clear that a variety of physico-chemical phenomena are associated with the 
polymerization of epoxy resins in an industrial reactor. These are quite complex and 
there is a definite need to study several of the interesting aspects of these 
polymerizations. Indeed, very few experimental or theoretical studies are available in 
the open literature on the modeling of industrial epoxy reactors . 5 " 7 This work intends 
to fill this gap to a small extent and attempts to provide a general mathematical 
framework for the study of the kinetics of the reactions. Limitations posed by heat and 
mass transfer effects can be incorporated subsequently, in order to develop a 
comprehensive reactor model. 

The first systematic study of the kinetics of epoxy polymerization was done by 
Batzer and Zahir . 1 ' 8 ' 10 These workers have reported both experimental and theoretical 
work on the molecular weight distribution of epoxy resins. Their work was based on 
Flory’s statistical method , 2 ’ 3 which, in turn, assumes equal reactivity of the functional 

Q Q 

groups. Deviations of experimental results from theory have been explained due to ’ 
the presence of monofunctional epoxides in the resin and the possibility of branching 
reactions. These workers conclude that the reactivity of epichlorohydrin with the 
phenolic hydroxyl group is higher than with the aryl glycidyl ether group. 

Ravindranath and Gandhi 4 improved the theoretical model developed by 
Batzer and Zahir and considered the unequal reactivity of epoxide groups. They used 
the probabilistic technique developed by Case" and observed improved agreement 
between model predictions and some of the experimental results of Batzer and Zahir. 
However, their model had two major drawbacks. It is well known that the alkali has a 


o 



very important role to play in determining the molecular weight distribution (MWD) 
of the product. Since alkali leads to the production of phenoxide ions , 1 it is possible to 
have some degree of control over the concentration of the phenoxide anion by 
controlling the rate of intermediate addition of the alkali. This enables one to alter the 
MWD of the resin. Unfortunately, the effect of alkali addition on the MWD is not 
considered by Ravindranath and Gandhi. In addition, in their model, one of the 
intermediate reactions ( dehydrochlorination ) has been considered to be extremely 
fast and hence the concentration of the chlorohydrin species is zero . This is contrary 
to the experimental observations of Batzer and Zahir , 8 who do observe this species in 
their GPC studies. Moreover, the use of the probabilistic approach by Ravindranath 
and Gandhi 4 prevents its use in modeling industrial epoxy reactors which are often 
operated in the semibaleh mode (with intermediate addition of sodium hydroxide and 
vaporization of the volatiles ). This study attempts to overcome these drawbacks, and 
presents a general kinetic framework, which can easily be adapted subsequently to 
incorporate several important physico-chemical phenomena which are often present in 
industrial reactors. 
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Chapter 2 


Formulation 

Any attempt to develop a mathematical model to describe the performance of an 
industrial semibatch epoxy reactor should first focus attention on the kinetics of the 
reaction, unencumbered by other influences. The detailed kinetic scheme 1,5 for epoxy 
polymerization incorporating most of the important reactions is shown in Table 1. The 
different end groups, -A,-B,-E and -F, are defined in this table. This scheme is more 
detailed than that used by Ravindranath and Gandhi 4 and takes care of the major 
weaknesses of their model, as described earlier. A study of the several reviews 5 ' 7 in the 
open literature suggests that the scheme in Table 1 is fairly general, and so should be a 
good starting point for the study of epoxy reactors. 

The various molecular species including the monomers present in the reaction 
mass at any lime, are shown in Table 2. AAq represents one of the monomers, bisphenol 
A, while EP represents the other monomer. This kind of description is required in order to 
have correspondence with the kinetic scheme of Table 1. Table 2 also gives the molecular 
masses of all these species ( for use in the computation of average molecular weights ). 

The kth moments (k=0,l,2,...) of the molecular species are defined in Table 3. It 
is necessary to use (n+l) k as the factor with the concentration, [i n ], of any species of 
‘length’ n to avoid giving zero weightage to the smallest species. The number and weight 
average molecular weights, M n and M w , for this system have been obtained starting from 
the fundamental equations 2,3 
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Table 1 

Reaction Scheme for Epoxy Polymerization 


1. 


OH + NaOH !— > -ONa + HoO 


A + NaOH — — — » -43 + W 


2. — ONa + CH 2 CHCH 2 C1 

V 

O 

- B + EP 


k 2 /H 2 0 


4 — OCH 2 CHCH 2 CI + NaOH 


k 2 /H 2 0 


OH 

F 


+ NaOH 


3. - 0CH 2 CHCH 2 C1 + NaOH — OCH 2 CHCH 2 + NaCl + H^O 

i V 

OH O 


+ NaOH 


k 3 


E + NaCl + H 2 0 


4. - ONa + CH 2 CHCH 2 0- 

V 

o 

-B + E- 


k 4 /H 2 0 


4 -OCH 2 CHCH 2 Q- + NaOH 
ill 


2 v_r 

l 


k,/H,0 


4 -(BE)- + NaOH 


5. — ONa+ C1CH 2 CHCH 2 0 ^-4 — 0CH 2 CHCH 2 O- + NaOH 

i ' 

* OH 

k 5 


— B + <-F 


OH 

-4 -<BF)— 


+ NaOH 


End Groups Present in the Reaction Mass 

'-A 

-OH 

- B 

-ONa 

-E 

— CH 2 qHCH 2 


0 

-F 

— OCH 2 CHOHCH 2 Cl 




Table 2 


Different Molecular(PoIymeric/Mononieric) Species Present in the 
Reaction Mass and Their Molecular Weights 

AA n : H(OROCH 2 CHCH 2 )„OROH 

OH 

AB„: H(OROCH 2 CHCH 2 ) n ORONa 
OH 


AE„ : H(OROCH 2 CHCH 2 ) n OROCH 2 CHCH, 

I V ' 

OH O 

AF n : H(OROCH,CHCH 2 ) n OROCH 2 CHCH 2 Cl 
OH OH 

BB n :Na(OROCH|CHCH 2 ) n ORONa 
OH 

BE n : Na(OROCH,CHCH 2 ) n OROCH 2 CHCH 2 

V 

OH O 

BF n : Na(OROCH,CHCH 2 ) n OROCH 2 j:HCH 2 Cl 
OH OH 


EE n : ,0^01 1CI I 2 (OROCH 2 Cl 1C1 I 2 ) n 0R0CH 2 CHpi 2 
O OH O 


EF n : v CH,CHCH 2 (OROCH 2 jCHCH 2 ) n OROCH 2 CHCH 2 Cl 

OH OH 


V 

o 


(Table 2 Continued ... a) 



Table 2 ( Continued. ..a ) 


FF„: ClCH 2 CHCH 2 (OROCH 2 £HCH 2 ) n OROCH 2 (j:HCH 2 Cl 


OH 


T 

OH 


OH 



where. 


CH 2 

Rt=-H 5 Cfi-C-C 6 H — 

CH, 

n=0, 1 , 2 ,... for all cases 


(Table 2 Continued... c) 



Table 2 (Continued ... c) 

Molecular Weights of Different Species 

Maaii = (h+1)(Maao+Mep - Mhci)- M p.p + Mhci 

Mabii = n(MAAo+ Mep - Mhci) + Mabo 

Maeo= (n+l)(MAAo+ Mep - Mhci) 

Mafii = (n+ 1 )(Maao + Mep - Mhci)+Mhci 

Mb Bn = h(Maao+ Mep - Mhci) + Mbbo 

MsEn = (ii+IXMaao+Mep - Mhci) +Mabo-Maao 

MaFn = h(Maao+Mep -Mhci )+ Mabo+Mep 

MEEn= (n+l)(MAAO+MEP - Mhci) +Mep - Mhci 

Mefh = (n+l)(MAA0+ Mep -Mhci)+Mep 

MFFn = (n+1 )(Maao+ Mep - Mhci) + Mep - Mhci 


where, 


n=0,l , 2 ,... for all cases 


Maao = molecular weight of bisphenolA 
Mep = molecular weight of epichlorohydrin 
Mhci= molecular weight of hydrochloric acid 






Table 3 

k"' Moments of Various Molecular Species 

= X(»+1)*[A4J 

w =0 

$AB = X(»+l)*MB fJ ] 

n=Q 

OO 

£ ae = Z(« + 1 )*M£„] 

/;=() 

^ af =Y J (n+\) k \AF n ] 

«=0 

A a /#/i = Z(// + !)*[££„] 

//=() 

^ he = Z<» + \) k \HE n \ 

//=() 

oo 

AV = £(n+l)*[fiFJ 
# 1=0 
OO 

**££ = Z<« + !)*[££„] 

# 1=0 

^ /'./■ = Z («+!)*!/?£„ ] 

n~() 

Xrr =Z(/i+l)*[FF„J 

n=0 


where, 

k=0, 1 , 2 ,... 


(Table 3 continued... a) 



Table 3 (Continued... a) 


M, 


Average Molecular Weights 

l ( M A'to + M /./> — M Hc i } { A 1 /\,a + A 1 Mt + A 1 ,\e + A 1 a i- + A 1 bb + A 1 ba + A 1 ba + 

A 1 at? + A'/r/.- + A 1 ff — A° a/? - A ° nn - A 0 ba } + Af £/ > {A°aa + A 0 /./.- + A° at-- + A°ba 
~A° A4 + [£P] 5 m0 } + M HCI {A 0 a t + A 0 44 + A°/\a - A°ae } + Af ABn {A 0 4« + A%£ 
_ +A°/;a } - M mo ?L rf. + A/ /;/} 0 A°/; fi ] 

{A [) AA + A () Alt + A°4A' + A°4£ + A° /{/{ + A%a + A%a + A°A't' + A 0 AT- + A°aa 

+[£/ , ] 5 ,„„) 


1 { A/ 440 + M FP - M ad } " { A - 44 + A' 4 n + A - 4 a + A" 4 a + A" bb + A ■ /;a + A" ba 

+A“ AA + A" AA' + A" FI- } + 2{ M AM) + M F p — M HCI } { ( Af hci — M HP ) A 1 44 + 

ABO + M HCI ~ A ^440 ~ Af aa)A 4 B + Af HCI AF + ( A fflBO + ^ HCI ~ ^ AAO 
-M AA )A III! + ( Af A/in - M AM) ) A BA +(M HP + M HCI + A/ A/)() - A/ mo - 
M hp ) A 1 nr + ( M E p ~ A/ hci ) A 1 aa + Af £/ > A 1 aa + ( A/ £/ , + M hci ) A 1 aa ) + 

(M HP - M HCI ) ~ ( A ° 44 + A° AA ) + ( Af 4 B0 + Af ffd — Af 440 - Af £/ > ) 2 A° 4 fl + 

A t af + ( Af B A 0 + Af HC/ - Af AA0 - Af £P ) 2 A°bb + (M abq - M AM )~ 

A 0 ba + ( Af aa- + Af HC , + Af Aeo - Af AA0 - Af EP ) A 0 ba + Af £ p 2 A 0 aa + ( Af EP 

+Af £C7 ) 2 A 0 AA +[£P]M 2 AA(5 m0 ] 

[{ Af AA0 + Af £p — M HCI } ( A 1 44 + A 1 4 B + A 1 4A + A 1 4A + A 1 Bfl + A 1 be + A 1 ba + 

A 1 aa + A 1 aa + A 1 aa — A° 4 /j — A° mi — A° b/.- } + Af /;/ , { A° aa + A° /./. + A° aa + A (1 ba 
-A° 44 + [ EP](5„, () } + Af HCI { A° aa + A°44 + A°4A - A°aa } + Af 4B0 { A°4B + A° ba 
A 0 ba } _ Af 44 0 A°ba + Af bb ( )A°bb 1 


PD/ = 


Af„ 


where, 

S m{) =0 for molecular weight determination without epichlorohydrin 
= 1 for molecular weight determination with epichlorohydrin 


in 



M n = ^(number fraction of i 111 species) x( molecular weight of i th species) (I) 

all species 

M w = ^(mass fraction of i 1 * 1 species) x (molecular weight of i ltl species) (2) 

all species 

The final expressions for these are also given in Table 3. Two sets of average molecular 
weights have been described. In one, the contribution of the epichlorohydrin is 
incorporated, while in the other, this is excluded from the computation. The latter is more 
appropriate and useful since experimental measurement of average molecular weights is 
usually performed after distilling out most of the epichlorohydrin from the reaction mass 
(there is negligible bisphenol A after some reaction has occured). 

The mole balance equations for a well mixed batch reactor for the various species 
shown in Table 2 can easily be written, and are shown in Table 4. In these equations, 
extreme caution has been taken to account for all the possible formation and consumption 
terms of the individual molecular species. The detailed reactions between the several 
molecular species required to develop the equations in Table 4, are given in the Appendix 
(Table A-l). The mass balance equations in Table 4 can be summed up analytically (using 
the definitions in Table 3) to give the ordinary differential equations (ODEs) for the 
several moments (k=0, 1 and 2). These are given in Table 5. The concentrations of the four 
end groups can be written in terms of the concentrations of the individual molecular 

species. These are given in Table 6. Mole balance equations for the various end groups 
and for some of the lower molecular weight oligomers (for which experimental data 1 are 
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Table 4 


Mass Balance Equations for The Molecular Species 


d[AA m ] 




dt 

d[ABJ 

dt 


i=0 


HI- 1 

f=0 


-'mO 


= *, [N«0//l{2[MJ-tAB m ]}-A : l/tBJ[EP]-^[/lBJ{£{[A£„] + [S£J+2[ EEJ + IEFJ} 

n =0 

-*,[ 4/t„, kXuafj + i /?/•;, i + 1 i + 2 | /•/•; i j + { 2 *., §i bb, n ab,„ , , 1 t k, Xl /w*. I! /?/•:,„ , , I 

n 0 i 0 »-(> 

//i ! /H I 

+ 2*, II AB, H BF m .,. x 1 + 2k, £l /«, ]( 1 


d[AE m ) 

dt 


= -k ] [AE„,][NaOH} + k,[AFJ[NaOH)-k 4 [AEJ&{[AB„] + [BE n ] + 2{BB n ] + [BF,,]} 


m~] 


m-\ 


n -0 

fit -l 


m • 1 


+ {ft4l[5£ 1 -][AE m .,.,J + 2fc 4 X[^B,][EE (I( . M ]4* s X[BE f J[AF B . ( . I ] + A ;i 2lAfi i I[£F OT . M ]}5 ll 


<=0 


i-Q 


d[AFJ 

dt 


= ] + *,[A5J[FF]- *,[AFJ[Mi<W]- M AFJ{£{ [AB, ] + [ MJ 


«= 0 
wi-I 


/=<) 


i-O 

+ 2* s £[AB,HFF M _ < _ 1 l}5 in0 




dt 


= k l [AB m }[Nc,OH}-2k,\BB m }\EP)-2k 4 [BB„,)[J J {[AE n ] + [BE„] + 2[EE„] + \EF„\} 


OO HI - 1 W 1 

2MFFJ{S{[^^,] + [^,l + 2[F^J + [^J} + {2/: 4 S[fiFJ[BE„,. | ..J + 2^X[/?EJl/iE„ 1 _,_a}<5, l , 


/t~0 


/=() 


/■=() 


^[BFJ 


, ^,[/lE ra ][^0//]-As[B£J[EP] + ^[Sf;j[/VflOH]-A: 4 [i?£ m ]{Xu^EJ + Ufi„] 

at 


;j=0 


+ 2| /i«„ 1 + 2| /i£ „ ] + 1 0FJ + 21 /■;/:„ ] + [ EF„ ]}-k,\ BE m ] { £ { l AF„ ) + \ BF„ \ + [ EF„ \ + 2( FF n ] ] 

n-Q 

m - | w 1 'M_l '»} 

+ {kiU BE, 11 BE,,, , , 1 + 4A: 4 Z( BB, ][ EE m _,_ , I + k< X[ BE, H BF„, , , | + 2*, 2,1 BB, H EF„, , , 1 )<5„ 


(Table 4 Continued ... a ) 
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d[BF m ] 


dt 


k 1 [AF m ][NaOH] + 2k 2 [BB m ][EP]-k 2 [BF m ][EP]-k 3 [BF m ][NaOH]-k4[BF m l{£fAE n ] 


n=0 


+ 


lBE] n +2[EE n ] + [EF n ]}-k 5 lBF m J{XllAF n J + lAB n J + lBE n j + 2lBB n ] + 2iBF n l + IEF n ] 


m-1 


n-0 


m-1 


m-1 


+ 2[FF n ]} + {k 4 £[BF,][BE m _ i _ ! ] + 2k 4 S[BBJ[EF m _ M ]+k 5 SfBFJtBF,^..,! 


i=0 


i=0 


i=() 


+ 4k 5 X[BB i ][FF m _ i _ 1 ]}8 m0 

i=0 


d[EE„ 

dt 


= k,[EFJ[NaOH]-2k 4 [EE m ]{X{[ABjnBEJ + 2[BB n ]HBFJ) + {2k 4 y Z[BE i ][EE m _ i _ { ) 


«=0 


nt-\ 


+ t s XlB£,][£F..,.,] !«,„ 


rfjEgJ 

dt 


= -k,[EFJ[NaOH] + k 2 [BEJ[EP) + 2k,[FF m ][NaOH]-k 4 [EFJ{'Z{[AB n \ + 2[BB n ) 


HBEJHBFJ}-k,[EFJ{2{[ABJ + 2[BBJ + {BEJ + [BF l] )} + {2k 4 J 4 [BF i }[EE m , i _ 

n=0 

in - 1 w-| »i- 1 

+ K Si BE, ][ EF m _ ] + k, Si BF, ][ EF m ._ t ] + 2k, Si BE, ][ FF m i ,})8„ 


i-0 


i-0 


d[FFJ 

dt 


= -2A-J FFJ[ NaOH) + k 2 [ EP}\BF m ] - - 2k, [ FFJ { £ 1 1 AB„ ] + 2[ BB n ] + 2[ BE n ) + \BF„}} 


n — 0 


MI-1 


+ (*, 2J Bf; ][ 1 + 2*, El BF, || FF. .,.,1)5. 


5 m o = 0, for m = 0 
5 m o = I , for m > 0 
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Table 5 


Moment Equations in a Batch Reactor 

dAp 

— ^ = -2k 1 [NaOH]A° AA + k 4 A? AB ^° AE + AB ^P AF 

dA? 

-^S- = -k 1 [NaOH]A. 0 AB+2k 1 A 0 AA [NaOH]-k 2 A 0 ABfEP]-k 4 A, 0 AB{^ 0 AE+2A°EE + >' 0 EF} 
- k 5 ^° AB {^°AF + EF + 2A° ff } + 2k 4 A° B B^° AE + 2k 5 A? BB^AF 

d/P a p 

-jf- = -*, [ NaOH]X° AE + k 3 X° Ap [NaOH ] - AE {X° AB + 2 A 0 BB + X° BF } + 2 k^ AB /9 EE 

+ k 5 $BE$ AF + k 5 $ AB^EF 

= -k x [ NaOH]X° AF + k 2 ft ab [EP] - k 3 ft AF [ NaOH) - A F i $ AB + ^ BB + $ BE } 

+ k 4 ^ BF$ AE + k 4$ AB $ EF + 2k 5 $ AB$ FF 

= k x {NaOH}$ AB - 2k 2 X° BB [EP)-2k 4 X 0 BB {X° AE +2X 0 EE + A 0 EE } - 2 k^ BB {A 0 AF 

+ A 0 EE +2$ EE } 
c/A° 

L -jp- = k { [NaOH)X 0 AE -k 2 $ BE [EP) + k 3 $ gplNaOm-k^BE^M + 2/9 AB + 2A° BjB 

+ $BE + $ BF + 2A° EE +A°/7/7}- k 5 i ] BE { A°af + $ EF + 2A° pp } + 4/^A 0 BB $ EE 

+ 2k 5 ?P BB /P EE 

= k { [NaOH]/9 AF + 2Jc 2 X° BB [EP)-k 2 ft BF [EP]-k 3 ft Bp [NaOH]- k^ BF {A° A £ 

+ 2/9 EE + A°£/r } - k 5 A 0 BF {$ AF + $ AB + 2A°f£ + $ BE + A°/j/r + A 0 EF + 2A° FF } 

+ 2/c 4 A°5bA° EF +4^ 5 A° BB p9 FF 

— = k 3 [NaOH] A 0 EF - 2k 4 A° EE { $ AB + 2A ° BB + /9 BF } + k^ A 0 BE A 0 Ep 

(Table 5 continued ... a) 
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Table 5 (continued ... a) 


dA o 

= k 2 [EP]A?BE - k 3 A° EF [NaOH] + 2k 3 A° FF [NaOH]- k 4 X 0 EF {^°AB + 2A () BB + ^BF} 

- k 5 X? ef( 2- 0 AB + 2 A? BB + A-°BE 1 + 2k 4 A° BF A° ee + 2k 5 A° B eA? FF 
= k 2 [EP]X°BF - 2k 3 X° FF [NaOH] - 2k 5 A° FF { A° AB + 2A° BB + A° BE } + k 4 X 0 BF X 0 E F 

= ^[NaOHjA 1 AA + k 4 {A ! AB A° AE + A 0 AB A ! AE } + k^A 1 AB A° AF + A 0 /^’ AF> 

= -kjtNaOHjA. 1 AB + 2k^ AA tNaOH] - k^ 1 AB [EP]~ k^ 1 AB {A° AE + 2A° E E + 

^EF ) ~ k 5 ^ AB {*? AF + A°EF + 2 ^° FF) + 2k 4 { A 1 BB^°AE + ^°BB^ AE) 

+ 2k 5 { A^bA 0 AF + A 0 bbA- 1 AF } + k 4 ^° AB^ BE + k 5 ^° AB^ BF 

^^• = -k 1 [NaOH]A , A E + k3X 1 AF {NaOH]-k 4 X 1 AE {A 0 A B+2X 0 BB+A 0 B F } + k 4 A. , BEA° A E 

+ 2k 4 {A. 1 AB A.°ee + A°A B A. , EE } + k 5 (A 1 beA°AF + A 0 beA 1 AF } + k 5 (a'abA 0 EF + A°abA ! EF 1 

= -k,[NaOH]A , AF + k 2 A ! AB [EP]- k 3 A ] AF [NaOH]- ^^{A 0 ^ +2A° BB + A° BE } 

+ k 4 (A 1 BF A° AE + A° BpA 1 AE ! + k 4 {A 1 AB A° E F + A f) AB A 1 FF } + k^A 1 BF A° AF 
+ 2k 5 { A 1 AB A° FF + A^gA 1 FF } 

kjtNaOHIA 1 ab - 2^ BB [EP] - 2^88 (A° A E + 2A° E E + A° EF ) " 2^' bb{A° AF + A° HF 
+ 2A° FF } + 2k 4 A°BBA , BE + 2k 5 A° BB^BF 

kjfNaOHIA 1 ae - ^A 1 BE [EP]+ k^ 1 Bp [NaOH]- ^a'be^ae + A° AB + 2A° BB + A° BB + A° BF 
+ 2A°ee + A 0 EF } - k 5 A , BE (A°A F + A° EF + 2A° FF } + k 4 A® BE A* BE + 4k 4 (A* BB A^ E E + A 0 bbA'ee) 
+ k 5 A° beA 1 bf + 2k 5 (A 1 BB A° EF + A° b B A ! EF 1 

= k^NaOHlA'Ap + 2k 2 A ! BB (EP]- k^ 1 Bp [ EP ] — k 3 >J Bp [NaOH]- k 4 A* B F { a 5 } A E +2A () E E 
+ 2.°EF) - k 5 A ] B F ( A°AF + ^°AB + 2A° BB + ^°BE + 2° BF + *°EF + 2A° FF } + k 4 A° BF A ! BE 
+ 2k 4 {A 1 BB A°ef + A°bb A 1 hi; } + k^BpA 1 B F + 4k 5 (A 1 BB^FF + A () B B^ FF 1 

(Table 5 continued ... b) 
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Table 5 ( Continued ... b) 


" ' d f E = k 3! NaOM IX 1 pp - 2k 4 ^ pp { X° AB + 2X° BB + X° BF ) + 2k 4 x' B p*-°pp * k 3 (X’ ppX^.p 
+ ^°BE^EF* 

= k 2 [EP]A?BE - k 3 x' FF [NaOH] + 2k 3 ^ ! Fp [NaOH] - k^pj^AB + 2X°B B + X°BP} 

- k^HFf^AB + 2X° BB + 2X° BE } + 2k 4 {X 1 BF X° EE + xVx'be} + + k^WW 

+ 2k 5 {X 1 BE^FE + X° BE^ 1 FF } 

= k 2 [EP]X 1 BF - 2k 3 X ! FF [NaOH] - 2k 5 X , FF {X° AB + 2X° BB + X° BE } + k 4 {X 1 BF X° EF + 
x° BF X* EF } + 2k 5 x' BF X° pp 

= -2k 1 [NaOH]X 2 AA + ^AsAe + A,iAe + AeAe ) + ^ AbAf 

+ X° AB X 2 AF + 2X 1 AB X* A p} 
dX 2 

— ^ = -k^NaOH^AB +2k 1 X 1 AA [NaOH]-k 2 X 2 AB [EP]- k 4 X 2 AB {X 0 AE + 2 X°be 
+ 2° EF } - k 5 X 2 AB {X° AF + X° ef + 2X° F F ) + 2k 4 {X 2 BB X° AB + X° BB X 2 AE 
+ 2X 1 bb^ae } + 2k 5 {X 2 B B^°AF + ^BB^AF + 2 X , B b2- 1 AF} 

+ k 4 {X° ab ^ 2 BE + 2 X 1 AB^ BE 1 + k 5 {^° AB ^ 2 BF + 2 X 1 A B^ BF 1 

^p- = -k 1 X 2 A E[NaOH] + k 3 X 2 AF [NaOH]-k 4 X 2 AE{ X° AB +2X° BB +X° BF } + k 4 {X 2 BE X° AE 

+ 2X 1 se^ab } + 2k 4 (X 2 AB X° EE + X° AB X 2 EE + 2X 1 AB x' EE } + k 5 (X 2 BE X° AF + X° BE X 2 AF 
+ 2X 1 bbX 1 AF } + k 5 {X 2 AB X° E p + X° AB X 2 EF + 2X 1 abX 1 EF } 

Pjp- = -k 1 [NaOH]X 2 AF + k 2 X 2 AB [HP]- k 3 X 2 A p|Nat)Hl- k 5 X 2 AF (X° AB + 2X° BB + X () BF .I 
+ k 4 {X 2 BF X° AE + X° BF X 2 AE + 2X 1 B pX' AE } + k 4 ( X 2 AB X° EF + X° AB X 2 EF + 2X 1 AB X 1 EF } 

+ k 5(X 2 BF X 0 A p + 2X 1 BF'X’aF I + 2k ? {X 2 AB X° FF + X° AB X 2 pp + 2X ! AB X , FF } 

^jp- = k,[NaOH]X 2 AB - 2k 2 X 2 gg[EP] - 2k 4 X 2 BB {X° AB + 2X° EE + X° EF ) - 2k 3 X 2 B B (X° AF + X° EF + 

2X° FF } + 2k 4 {X°B B X 2 BE + 2X 1 BB^BE ) + 2k 5 {X° BB X 2 BF + 2X 1 BB X 1 BF ) 

(Table 5 continued ...c) 
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Table 5 ( Continued ... c) 


cR 2 

— ^ = k,fNaOH]A. 2 AE - k 2 A. 2 BF i Ep J + k 3^ 2 BF I NaOH] - k 4 A. 2 BH { A.° AR + X° AB + 2 X° BB + A°be + X% 

+ 2A.° ee + A°efI ~ k 5^ 2 BE f^°AF + ^°EF + 2A.° FF I + k 4 {A.°BE^ 2 BE + 2A. 1 BE*- 1 BE ! 

+ 4k 4 {A. 2 bb^ 0 EE + A? BB^ 2 EE + BB^ EE } + k 5 { A°beA- 2 BF + 2A. 1 be*- 1 BF 1 
+ 2k 5 {X? bb^° EF + A? BB^EF + 2 A. 1 BB^’ EFl 

= k ] [NaOH]A? AF + 2k 2 X 2 BB (EP]- k 2 A? BF [EPl- k 3 X 2 Bp [NaOH]- k 4 X 2 BF {X° A E + ^EE 

+ ^’hf) - k 5^- 2 BF^AF + A°ab + 2X°bb + ^°BE + BF + A. U EF + 2 ^°FF 1 + k 4 {^°BF^ 2 BE 
+ 2A-*BF^BE} + 2k 4 {A. 2 bb^- 0 EF + 2.°BB^ 2 EF + 2 ^ BB'’ 1 ' EF) + k 5^BF^ 2 BF + 2 ^BF^BF) 

+ 4kg {A 2 BB^°FF + ^°BB^ 2 FF + 2 ^' BB^'ff) 

2 

= k 3 lNaOH]X 2 EF-2k 4 A. 2 EEf^AB + 2 ^°BB + ^Bp) +2k 4 {2- 2 BE^EE 4-2 ^ BE^EE* 

+ k 5 BE EF + ^° BE ^ EF +2 ^ BE ^ EF ^ 

dA 2 

-^f E = k 2 lEP]A. 2 BE " k 3 A. 2 EF [NaOH] + 2k 3 X 2 FF [NaOH]- k 4 X 2 EF{X°AB + 2A? BB + A-Vl 

- k g A 2 EF {A? AB +2 X°bb + 2aP B E} + 2k 4 (A 2 BF A.° EE + A.° BF A 2 EE + 2A. 1 BF A. 1 ee} + k 4 (A. 2 beA°HF 
+ 2A. 1 be^ 1 EFl + k 5 (A 2 BfA-°EF + 2 A. 1 BF A.' EF } + 2kg {A 2 B eA-° FF + A- 0 BE A 2 FF + 2A 1 BE aJ FF’} 

= k 2 [EP]A. 2 BF - 2k 3 X 2 FF [NaOH) - 2kgA 2 FF f A° AB + 2A° BB + A. () BB } + k^fiF^EF 
+ A? BF A 2 ef + AA.'bfA 1 ef ) + 2kg (A 2 BF A.° FF + 2/J b F A.' FF } 
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Fable 6 


Concentration of Various End Groups 


[A] = X( 2[ AA n ] + [ AB n ] + [ AE n ] + [ AF n ] ) 
n=() 

on 

iB]= X(2[BB n ] + [AB n ] + [BE n ] + [BF„J) 

n=0 

oo 

[El = E(2[EE n ]+[AE n ]+tBE n ] + [BF n ]) 
n=0 

OQ 

(F|= Z(2[FF„]+lAF„]+[Bn n ] + [EF n ]) 

n=() 
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available) can be written using the simplified kinetic scheme (in terms of end groups ) of 
Table 1. These are given in Table 7. The following checks were made to confirm the 
correctness of the balance equations. 

(a) It was found that the equations obtained by appropriately summing up the 
equations in Table 4 (using Table 6), were identical to the balance equations for the end 
groups (given in Table 7 and derived independently using the simplified kinetic scheme 
of Table 1). The same check was made using the equations for the zeroth moments ( 
Table 5) and the equations matched with those in Table 7. This confirmed the correctness 
of the equations for the zeroth moment, as well as partially confirmed the correctness of 
the mole balance equations in Table 4. 

(b) In a batch reactor, the total moles of the aryl unit, R, is conserved, i.e., the moles 
of R at any time ( distributed over the several molecular species given in Table 2) must be 
equal to the moles of R initially (in bisphenol A). This has been checked out, and 
confirms that the equations are correct. 

(c) In a batch reactor, the sum of the moles of the two groups, — CHCH 2 and -CH 2 CH — 

O OH 

should be constant with time and should be equal to the initial moles of epichlorohydrin 
present. An equation for the total number of the above repeat units (distributed over the 
several molecular species in Table 2 ) was written ( in terms of the moments ). It was 
found that this check was also satisfied. 
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Table 7 


Mole Balance liquations for the bind Groups and Some Lower Oligomers 

in a Batch Reactor 

d[A] 

— — = -k\[NaOH\[A] 


d[B] 

dt 


k { [NaOH][A] - k 2 [EP)[B } - k 4 [E][B } - * 5 [F][fl] 


d[E] 

dt 


k 3 [NaOH][F]-k 4 [B\[E] 


d[F] 

dt 


k 2 [EP}[ B] - fc,[ NaOH][ F ] - Jfc 5 [ B\[ F ] 


d[EP] 

dt 


-k 2 [EP][B] 


= -k^NaOHllA] + k 2 [EP][B] - k 3 [NaOH\\F] + k 4 [B\[E\ = 

~ - y C -- = k,[NaOH][F] + k,[B\[F] 
at 

= -Ik^NaOH^AA^ 

= -k { [NaOH][ AE 0 ] + k 2 [AF 0 ][NaOH] - k 4 [AE 0 ][B] 


(Table 7 Continued ... a) 



Table 7 (continued ... a) 



--^[^O/fj[ylFj + *jAB 0 J[£P]-^[^ 0 ][^//]-^[AF 0 ][5] 



k 2 [EP][BF f) \ - 2k 3 [FF 0 ][NaOH] - 2k,[FF 0 ][B] 


d[EE 0 ] 

~dT~ = k * [ NaOH i EF o] ~ 2k 4 [EE 0 ][B] 


4ab 0 ] 


dt - -k,[NaOIi\ {[AZ? 0 ] - 2[AAJ } - k 2 [AB 0 )[EP] - k 4 [AB 0 j(£j - * 5 f AZ? () ]fFJ 


4 B ^o] 


(It k ' ^ N “OH][ aE () ] k 2 \BE 0 ][ EP\ + kj BF () ][ NaOH ] - k A [ BE () ]{ [ B \ + [ E] } - *,| BF t) )( F | 


d[EI^] 


dt 


— = *, 1 NaOH || AB () ] - 2 k 2 [ BB {) If EP ] - 2k , [ BB {) If E\ - 2k, f BB () |[ F | 
d[BFf)] 

= ki[NaOH][AFb] + 2k 2 [BB 0 ][EP] - k 2 [Bfb][EP] - k.^fB^JfNaOH] 

- k 4 [BF^)] { [ B] + [E] } — k^fBfylfFJ 

= k 1 [EP][BE 0 ]-k3[EFb][NaOH3 + 2k3[FF 0 ][NaOH]-k 4 [EF () ][B]-k 5 fEFf ) ]rB] 


-fc^TRAL LIBRARY 

I I. T.. KANPUR 


The above three checks were sufficient to confirm the correctness of the mole and moment 
balance equations in Tables 4 and 5. 

The equations in Table 5 ( alongwith those for [NaOHJ and (EP| from Table 7 ) 
form a complete set of ODEs (initial value problem-IVP), which can be integrated for any 
given set of initial conditions for a specified temperature history and specified rate 
constants (ki-k<0. Similarly, the equations in Table 7 form another complete set of ODE- 
IVPs and can be integrated separately. The numerical technique used to integrate these 
equations is Gear’s algorithm. 12 The NAG library routine, D02EBF, which has a built-in 
step size control algorithm and is particularly useful for stiff systems, is used (with a 
tolerance, TOL, of 10' 7 ) for this purpose. Changing the value of the parameter, TOL, led 
to almost identical results. The three checks described above were also made with the 
computer code to ensure the correctness of the simulation program. 

With the mathematical model and computer code now developed, we can start 
obtaining some useful results. Some experimental data on isothermal ( 25" C) 
polymerization of bisphenol A with epichlorohydrin with NaOH catalyst (in cellosolve 
solvent ) has been provided by Batzer and Zahir. 1 These could be used to obtain curve-fit 
values for the five rate constants at 25° C. Similarly, curve fitting can be attempted at other 
temperatures when additional experimental data becomes available to give Arrhenius 
expressions for the rate constants ( unfortunately, experimental work along the lines 
required are not available in the open literature ). The Box complex procedure 1314 is used 
to carry out the curve-fitting exercise for the data presented by Batzer and Zahir. 1 This 
technique minimizes an objective function, E, using a pattern-search technique. The 
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objective function is chosen as a weighted sum of absolute errors between the values 
predicted by the model and experimental data. It is given by 


E( k t ,k 2 ,k^,k 4 ,k s ) 


j 1N J,exp i 


[j] i,thcor -UJ i,exp 


[jli. 


,exp 


(3) 


In the above equation, [j] is the concentration of the jth molecular species for which 
experimental data is available, 1 Wj is the weightage factor assigned to the experimental 
data set for the jth oligomer (or summation of oligomers), and subscript i indicates the ith 
data point in a set for the jth species. Superscripts, exp and theor, represent experimental 
and the corresponding theoretical values for any concentration. Nj. CX p indicates the number 
of experimental data points available for the jth species. Minimization of E while 
searching for values of kj-ks leads to a nonlinear best-fit set of values for these rate 
constants. The CPU time taken on a mainframe HP9000/850 S supermini computer for the 
parameter estimation (220 iterations) is 28.64 s. 
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Chapter 3 

Results and Discussion 

The Box complex algorithm is used with the experimental data of Batzer and 
Zahir 1 to obtain the best-fit values of the five parameters, k,, k 2 , k 3 , k 4 and k 5 , at 25° C. 
The initial conditions used are : [AA 0 ] 0 =200 mol/m 3 , [EP] 0 =1280 mol/m 3 , [NaOH] 0 =400 
mol/m 3 . The parameters used in this program are given in Table 8. The experimental 
data' used in our study are the concentrations, [AA 0 ] (bisphenol A), [EE 0 ], [EF 0 ], [FF 0 ] 
and ([AEo]+[AFo]), at different times. The actual data of Batzer and Zahir ( representing 
uncorrected peak heights of chromatographic curves ) were nondimensionalised using the 
corresponding maximum values, while the chromatographic peak data on [AA 0 ] was 
transformed into concentrations using the given initial concentration, [AA 0 ]o- The optimal 
values of the rate constants are given in Table 8. 

Fig. 1 shows how E varies with iteration number. It is observed that E slowly 
reaches its lowest value. Figs. 2, 3 and 1 show how significant changes occur in the 
values of kj-k 5 before they also settle down to their best-fit values. 

The best-fit values of the rate constants are used to generate the concentration 
histories of several lower oligomers. Figs. 4 - 6 show the quality of the curve-fit of the 
experimental data. These fits are observed to be fairly good. Figs. 7 and 8 show the 
variations of [EP], [NaOH], [NaCl] and [H 2 OJ with time, computed using the best-fit 
values of the rate constants. It is observed that reasonable amounts of NaOH are 
consumed leading to the production of water. The epichlorohydrin concentration does not 



Table 8 


Data Used and Results Obtained in Optimal Parameter Estimation 
Weightage Factors Used in Eq. 3 


Data Used 

Wj 

[AA„] 

1.0 

[EE 0 ]/[EE 0 ] max 

1.0 

[EF 0 ]/ [EF 0 ] max 

1.0 

[FFo]/[FF 0 ] max 

5.0 

([AE 0 ]+[AFo])/([AEo]+[AF„]) max 

1.0 


Values of the Parameters Used in Box Complex Optimization 

N=7 
a =1.3 
(5 =1()- 6 
y=3 

5 =0.01(hj-lj) where, hj and lj are the upper and lower bounds on the jth parameter 
Initial Value of the Parameters 

[k, ,k 2 ,k 3 ,k 4 ,k 5 ] = [1.2x 10 -3 ,3.1x 10~ 5 ,2.0x 10" 4 ,6.4x 10 _6 ,2.7x 10 1 ] 


Bounds on the Parameters 

5.0 x 10 -4 <k, < 15x 10 -3 
l.OxlO -5 <k 2 < 9.0 x 10 -5 
l.OxlO -4 <k 3 <9.0xl0 -4 
l.OxlO -6 <k 4 < 9.0 x 10~ 6 
l.OxlO 1 <k 5 <9.0 x10 ! 


Optimized Values of the Parameters 


kl 

k2 

k 3 

k 4 

L k 5 J 


1.01 1 378 xlO -3 
6.31 6743 x 10 -5 
3.949022 XlO -4 
2.232 1 728 xlO" 6 
3.31 16070x10' 




Fig. 1 : Variation of error and k 5 with iteration number. 
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Fig. 4 : Variation of the concentration of bisphenol A with time 
Experimental data of Batzer and Zahir 1 also shown. Curve 
indicates model predictions using the optimal parameters. 



0.8 



time, hr 

Variation of the dimensionless concentrations of EE 0 andFF 0 
with time. Solid lines represent model predictions using optimal 
parameters. Symbols represent experimental data points of 
Batzer and Zahir 1 . 
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Fig. 7: Variation of [EP] and [NaCl] with time. Optimal values of the 
rate constants used. 





change much with time since it is taken in excess. Small amounts of NaCl are produced 
during the polymerization. This needs to be removed by washing with water at the end. 

The optimized model parameters are used in the moment equations of Table 5 
with [AAo]q— 200 mol/m' , |EP]o=l 280 mol/m 3 and [NaOH]o=400 mol/m 3 , to generate the 
n umbei and weight average molecular weights, as well as the polydispersity index (with 
the help of Table 3 ). Fig. 9 shows the variation with time of the average molecular 
weights with and without the concentration of epichlorohydrin used in the computations. 
The value of M n (without EP) is slightly over 342, the molecular weight of DGEBPA, 
near the end, and indicates that the end-capping of bisphenol A by excess epichlorohydrin 
is essentially complete. The PDI without EP rises from 1.0 to about 1.4 during this 
period. T he PDI of about 1.4 at the end suggests that significant amounts of higher 
oligomers arc also formed in this process. 

Detailed sensitivity tests are now performed to identify the most important 
parameters, so that greater care can be taken in estimating their values. Fig. 10 shows that 
the bisphenol A concentration is somewhat sensitive to changes in the value of kj. The 
other rate constants do not influence the concentration history of [AA] {) as much. Similar 
sensitivities are observed for the concentration of NaOH. The concentration of NaCl, in 
contrast, is sensitive to variations in both ki and k 2 , but not as sensitive to other rate 
constants. The epichlorohydrin concentration is not influenced much by any of the rate 
constants, since it is in large excess. The value of M n ( without EP ) is most sensitive to 
variations in kj and k 2 . Fig. 1 1 shows how higher values of M n are attained when kj or k 2 



Variation of M n and M w with time. Dotted lines indicate average 
molecular weights incorporating epichlorohydrin. Solid lines 
indicate the variations excluding epichlorohydrin in the computation. 
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IS increased by about 20% ( over the optimal or reference values ). The sensitivity of M n 
to the remaining three rate constants is much less. Similar conclusions are obtained for 
the sensitivities of M w and PD1 to variations in the rate constants ( detailed plots can be 
provided on request ). A word of caution may be mentioned at this stage. The kinetic 
scheme of I able 1 indicates that the main chain-lengthening steps are reactions 4 and 5. 
The telative insensitivity of M n and M w to variations in k 4 and ks could possibly be 
because of enoneous optimal values of these two rate constants (experimental data on M„ 
or M w ‘lie not available, and the curve-fit using data on only the low molecular weight 
molecular species could bias our optimization against good values of k 4 and k 5 ). 

An interesting simulation study is now made. Several industrial operations 
involve intermediate addition of NaOlI. We run our simulation code to sec the effect of 
sudden changes in NaOH concentration. The results so generated would give us some 
idea about the effects of NaOH addition in industrial epoxy reactors. Fig. 12 shows that 
higher M n (and higher M w and PD1 ) are obtained with intermediate addition of NaOH. A 
similar effect would be expected if we distill off water/epichlorohydrin in between ( the 
concentration of NaOH would increase ). It may be added that the probabilistic model of 
Ravindranath and Gandhi 4 cannot be used to provide such results and that our kinetic 
model is quite general. If we use appropriate values of the rate constants ( corresponding 
to higher temperatures ), and initial conditions typically encountered in the taffy or 
advancement processes, we could easily obtain the results of Ravindranath and Gandhi. 
In fact, the latter model is really a subset of the present one. The probabilistic model. 
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Fig. 12 : Solid curve indicates the variation of M„ with time, with [NaOH] 
concentration increased by 400 mol/m 3 above actual concentration 
at 3 hours. Dotted curve indicates the variation without any 
intermediate addition of NaOH. 
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Chapter 4 


Conclusions 

A detailed kinetic model has been developed for epoxy polymerizations. The rate 
constants are curve-fitted using one set of experimental data available in the open 
literature. A sensitivity study has been carried out. The model incorporates the 
probabilistic model of Ravindranath and Gandhi 4 as a subset, but is far more powerful, as 
demonstrated by predictions for the case of intermediate additions of NaOH. 



Chapter 5 

Suggestions for Future Work 

Experiments should be carried out at different temperatures to study the 
temperature effect on the model parameters and develop Arrhenius equations for the rate 
constants. Experimental data are required on concentrations of the different oligomeric 
species as well as on average molecular weights, so that the model parameters do not get 
biased by the former. Thereafter, the kinetic model developed here can be extended to 
apply for industrial reactors, and can be tested (or tuned) against industrial data. 
Optimization studies can then be performed. 
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APPENDIX 


Table A-l 

Different Reactions Taking Place 

AA„ + NaOH AB n + H 2 0 

AB„ + NaOH — ^ BB n + H 2 0 
AE„ + NaOH — ^-4 BE,, + H 2 () 

AF„ + NaOH — > /?/;, + ///> 

AB„ + EP V " ;0 > AF; + EaOE 

BB„ + EP > BF;, + /ME 

BE,, + EP V " : ° > EF;, + AME 

BF; + EP - ° -> FF n + NaOH 

AF n + NaOH — > AE„ + AM + E 2 0 
BF;, + NaOH — > BE „ + AM + E,0 
EF;, + NaOH — > £E„ + AM + H 2 0 
FF n + NaOH - — > EF„ + AM + E 2 0 

AE„ + AE,,,-^^ AA + M 

BE,, + a£ ((i _H AB n+m+l + NaOH 

BE,, + AE„, . V*’ 0 - » A£, l+m+l + NaOH 

BF„ + AE,„ AF n+m+1 + NaOH 

AB n + BE m W > AB„ +m+1 + /ME 
BB„ + BE„, > BB„ + „ I+ , + /ME 

BE,, + BE,,, 2tj ///; — -» BE, l+m+l + /ME 
BF„ + BE. BF n+m+] + NaOH 


(Tabic A-l Continued... a) 



Table A-l : (Continued... a) 


AB n + EE m -*'"*° > AE n+m+l + NaOH 
BB„ + EE m > BE H+m+l + NaOH 

BE,, + EE m -* ,H '° ) £E„,„ (1 + NaOH 
BF n + EE m -! h ’"'-° > EF, , + NaOH 


AB n + EF m --*" h ° ) AF n+m+] + NaOH 
BB„ + £F,- 3 ' 3 ° > gF. w , + NaOH 
BE,, + EF m > £C + „, +1 + NflOH 
BF„ + £F m — ) ££_, + MzOff 


/2+/;h-1 


AB„ + AF m - 

*5 

AA,, +m+] + NaCl 

BB,, + AF m - 

2*5 

— > AB n+m+] + NaCl 

BE n + AF„- 

*5 

AE „ +m+ , + NaC/ 

BF+AF- 


-4 /\E + A/aC/ 

ti in 


iHHHl 

AB„ + ££„,— 

*5 

■> AB n+m+] + AtaCZ 

BB„+BF„- 


BB n+m+] + NaCl 

BE„ + BF„- 

*5 

*BE n+m+ , + NaCl 

BF„ + BF m — 

*5 

* B F n+m+] + NaCl 

AB n + EF m - 

*5 

-4 A£„ +m+1 + NaCl 

BB„+EF m - 

2*5 

-4 BE n+m+] + NaCl 

BE„ + ££„,- 

*5 

EE n+m+ , + AfaC7 

BF n + EF m - 

*5 

~>EF n+m+ 1 + AflC/ 


A£„ + ££ m -^Af;, +m+1 + yvaa 


££„ + ££„, 
BE„ + ££„, 
BF„ + ££„, 


4 *. 


*-» S' 7 ...., + WflC/ 


2 A- 


£F.„„ + NaCl 


2k, 


EF, + NaCl 


where 

m,n=0, 1 , 2,... 




